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Holographic vector mesons from spectral functions 
at finite baryon or isospin density 
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We consider gauge/gravity duality witii flavor for the finite-temperature field theory dual of the AdS- 
Schwarzschild black hole background with embedded D7-brane probes. In particular, we investigate spectral 
functions at finite baryon density in the black hole phase. We determine the resonance frequencies correspond- 
ing to meson-mass peaks as function of the quark mass over temperature ratio. We find that these frequencies 
have a minimum for a finite value of the quark mass. If the quotient of quark mass and temperature is increased 
further, the peaks move to larger frequencies. At the same time the peaks narrow, in agreement with the for- 
mation of nearly stable vector meson states which exactly reproduce the meson mass spectrum found at zero 
temperature. We also calculate the diffusion coefficient, which has finite value for all quark mass to temperature 
ratios, and exhibits a first-order phase transition. Finally we consider an isospin chemical potential and find that 
the spectral functions display a resonance peak splitting, similar to the isospin meson mass splitting observed in 
effective QCD models. 
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Recently in the context of gauge/gravity duality, there has 
been an intensive study of the phase diagram of J\f ~ 4 large 
Nc SU{Nc) supersymmetric Yang-Mills theory with added 
fundamental degrees of freedom, by considering the AdS- 
Schwarzschild black hole background with added D7-brane 
probes Ol H H H H H] . There are two kinds of D7-brane 
probes in the black hole background; Either they end before 
reaching the black hole horizon, since the wrapped by the 
D7-brane probe shrinks to zero as in JtII, or they reach all the 
way to the black hole horizon. The first class of embeddings is 
usually called 'Minkowski embeddings', while the second is 
referred to as 'black hole embeddings'. The parameter which 
parametrizes different embeddings is the temperature normal- 
ized quark mass niq/T, which may be given in terms of the 
asymptotic value xo of the embedding coordinate at the AdS 
horizon. The phase transition between both classes of embed- 
dings is of first order. The analysis of the meson spectrum 
shows that this phase transition corresponds to a fundamen- 
tal confinement/deconfinement transition at which the mesons 
melt. 

Particular interest has arisen in the more involved struc- 
ture of the phase diagram when a baryon chemical potential 
is present 181]. It was argued that for non-vanishing baryon 
density, there are no embeddings of Minkowski type, and all 
embeddings reach the black hole horizon. This is due to the 
fact that a finite baryon density generates strings in the dual 
supergravity picture which pull the brane towards the black 
hole. A chemical potential for these baryons corresponds to 
a vev Aq for the time component of the gauge field on the 
brane. In the dual thermal SU {Nc)-gauge theory a baryon is 
composed of Nc quarks, such that the baryon density n b can 
be directly translated into a quark density Ug = Ub Nc- The 
thermodynamic dual quantity of the quark density is the quark 
chemical potential jiq. In the brane setup we use, the chemical 
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FIG. 1: The phase diagram for quarks: The quark chemical po- 
tential yLq divided by the quark mass is plotted versus the tempera- 
ture T divided by M — 2mq/\f\. Two different regions are dis- 
played: The shaded region with vanishing baryon density and the 
region above the transition line with finite baryon density, in which 
we work here. The multivalued region at the lower tip of the transi- 
tion line is not resolved here. The curves are lines of equal baryon 
density parametrized by d = / {Nj^fXT^). The critical den- 

sity d* = 0.00315, at which the first order phase transition between 
two black hole phases disappears, is shown as short-dashed line close 
to the transition line. It virtually coincides with the short-dashed line 
for d = 0.002. 

potential is determined by the choice of quark density and by 
the embedding parameter xo ■ 

Very recently, however, it was found that for a vanishing 
baryon number density, there may indeed be Minkowski em- 
beddings if a constant vev is present, \yhich does not de- 
pend on the holographic coordinate 13, [ifl [lH [iSll- The 
phase diagram found there is sketched in figure [T| In the grey 
shaded region, the baryon density vanishes {ns = 0) but tem- 
perature, quark mass and chemical potential can be nonzero. 
This low temperature region only supports Minkowski em- 
beddings with the brane ending before reaching the horizon. 
In contrast, the unshaded region supports black hole embed- 
dings with the branes ending on the black hole horizon. In this 
regime the baryon density does not vanish {ub ^ 0). In this 
paper we exclusively explore the latter region. At the lower tip 
of the line separating ns = from ns 7^ in figure [T] there 
exists also a small region of multivalued embeddings, which 
are thermodynamically unstable ifioll . 

In the black hole phase considered here, there is a funda- 
mental phase transition between different black hole embed- 
dings H. This is a first order transition, which occurs in a 
region of the phase diagram close to the separation line be- 
tween the two regions with vanishing (gray shaded) and non- 
vanishing (unshaded) baryon density. This transition disap- 
pears above a critical value for the baryon density ns given 
by 

J* = 0.00315, d = 2^/^nB/{Nfy/XT^). (1.1) 

In this paper we make use of the methods developed in the 
context of AdS/CFT applied to hydrodynamics, for instance 
ifTii [Tsi [l6ll . in order to determine the spectral function at 
finite temperature and finite baryon density. For vanishing 



chemical potential, a similar analysis of the spectral functions 
has been performed in UtIi . It was found that the spectrum 
is discrete at large quark mass, or equivalently at low temper- 
ature. At low quark mass, a quasiparticle structure is seen 
which displays the broadening decay width of the mesons. 
As the mass decreases or temperature rises, the mesons are 
rendered unstable as the resonance frequencies develop imag- 
inary parts. Modes corresponding to such frequencies are 
called quasinormal. These excitations are then dissipated in 
the plasma. - Note that for this case, there are also lattice 
gauge theory results ifisll . 

In this paper we study the differences in the spectral func- 
tions with and without chemical potential. Relating our work 
to the phase diagram shown in [10, figure 2] (and reproduced 
here in figure [TJ, we here consider the region of black hole 
embeddings (unshaded region) with nonvanishing quark den- 
sity. We find that at low temperature to quark mass ratio, 
the spectrum is asymptotically discrete and coincides with the 
zero-temperature supersymmetric meson mass formula found 
in ifigll . which in our coordinates reads 

M =^y'2in+l){n + 2). (1.2) 

In lfl9ll . Loo denotes the asymptotic separation of the D3- and 
D7 -branes and n counts the nodes of the embedding fluctua- 
tions. Here we are considering s-wave modes in the Kaluza- 
Klein expansion of the D7 brane probe wrapping S'^, so the 
angular quantum number / is zero. We connect the structure of 
the spectra found to the phase diagram in figure[T] The meson 
mass behavior described above occurs close to the Minkowski 
region of the phase diagram, where temperature effects are 
subdominant. Moreover, as a function of decreasing temper- 
ature to quark mass ratio, the quasiparticle peaks behave dif- 
ferently with and without finite quark density. As discussed 
in [17], at vanishing density 71b ~ the peak maxima move 
towards smaller frequencies as a function of increasing quark 
mass. Here, in the case of finite quark density, we observe 
a similar behavior at small quark mass. However, keeping 
the temperature fixed as we increase the quark mass further, 
the peaks turn around at a value m'"" and move to larger and 
larger frequencies as the associated mesons become more sta- 
ble. Note that a turning point behavior was also observed for 
vanishing quark density in the context of quasinormal modes 
for scalar modes of melting mesons [@] . 

Our spectra also show that for given quark mass and tem- 
perature, lower n meson excitations can be nearly stable in the 
plasma, while higher n excitations remain unstable. At van- 
ishing baryon density, the formation of resonance peaks for 
higher excitations has also been observed in ||20|] . We discuss 
the distinct behavior of resonance peaks in section llllB 21 in- 
cluding a comparison of the observed turning points at finite 
baryon density with previous results. 

We also calculate the quark diffusion constant D and show 
that at finite density, it exhibits the first-order fundamental 
phase transition up to the critical density given by d* = 
0.00315. For very large values of the density, the diffusion 
constant asymptotes to D ■ T = l/(27r). This reflects the fact 
that in this case, the free quarks outnumber the quarks bound 
in mesons. 



3 



As a second point we consider the case of an isospin chem- 
ical potential, on which previous work in the holographic con- 
text has appeared in ll3[ l2lll . In this case, two coincident D7- 
brane probes are considered. In particular we extend the re- 
sults of our previous paper Ii22i1 . in which we calculated the re- 
tarded Green function and diffusion coefficient at finite SU (2) 
isospin chemical potential for the flat embedding uig = 0. In 
this previous work we also restricted to the case of constant 
vev for the non-Abelian gauge field Aq, where 3 is the flavor 
and the Lorentz index. This means that we chose to be 
independent of the AdS radial direction. In this case we found 
a non-analytic frequency dependence of the Green functions 
and the diffusion coefficient. Here we extend this work to the 
case of non-vanishing quark mass, leading to non-trivial D7 
embeddings, and to the case of radially varying gauge field 
component Aq. We find that spectral functions quantitatively 
deviate from the baryonic background case. Additionally, a 
splitting of quasi-particle resonances is observed, which de- 
pends on the magnitude of the chemical potential. 

This paper is organized as follows. In the following sec- 
tion [III we introduce the gravity background, field and brane 
configuration, used for the subsequent calculations. We also 
sketch the method to obtain retarded real-time correlators of 
thermal field theories from supergravity calculations. In sec- 
tion|III]we discuss the spectral functions and diffusion behav- 
ior of fundamental matter at finite baryon density. For matter 
with isospin chemical potential, the same analysis is carried 
out in section HV] The results are briefly summarized in sec- 
tion[Vl 



II. HOLOGRAPHIC SETUP AND THERMODYNAMICS 

A. Background and brane configuration 

We consider asymptotically AdS^ x space-time which 
arises as the near horizon limit of a stack of Nc coincident 
D3-branes. More precisely, our background is an AdS black 
hole, which is the geometry dual to a field theory at finite tem- 
perature (see e.g. 112311 ). We make use of the coordinates of ^ 
to write this background in Minkowski signature as 



with the metric dilg of the unit 5-sphere, where 

f{9)-i-4, m = i + ^\ ,,,, 

g Q (2.2) 

i?^ = 47r5siVca'^ qh = Tt:R^. 

Here R is the AdS radius, gs is the string coupling constant, 
T the temperature, Nc the number of colors. In the following 
some equations may be written more conviniently in terms of 
the dimensionless radial coordinate p ~ g/ Qh, which covers a 



range from p = 1 at the event horizon to p oo, representing 
the boundary of AdS space. 

Into this ten-dimensional space-time we embed Nf coin- 
ciding D7-branes, hosting flavor gauge fields A^^. The em- 
bedding we choose lets the D7-branes extend in all directions 
of AdS space and, in the limit p ^ oo, wraps an S^ on the 
S^. It is convenient to write the D7-brane action in coordi- 
nates where 



2 + sin2 9dnl), 
(2.3) 

with < 9 < tt/2. From the viewpoint of ten dimen- 
sional Cartesian AdS^ x S^, 9 is the angle between the sub- 
space spanned by the 4,5,6,7-directions, into which the D7- 
branes extend perpendicular to the D3-branes, and the sub- 
space spanned by the 8,9-directions, which are transverse to 
all branes. 

Due to the symmetries of this background, the embeddings 
depend only on the radial coordinate p. Defining x = cos 9, 
the embeddings of the D7-branes are parametrized by the 
functions x{p)- They describe the location of the D7-branes 
in 8, 9-directions. Due to our choice of the gauge field fluc- 
tuations in the next subsection, the remaining three-sphere in 
this metric will not play a prominent role. 

The metric induced on the D7-brane probe is then given by 



ds' 



R 



f 

,2 



!dx' 



1 



R\\-X^)dVLl 



■dp" 



(2.4) 



Here and in what follows we use a prime to denote a derivative 
with respect to q (resp. to p in dimensionless equations). The 
symbol \J~g denotes the square root of the determinant of the 
induced metric on the D7-brane, which is given by 



(2.5) 



The table below gives an overview of the indices we use to 
refer to certain directions and subspaces. 



(2.1) 


coord, names 


■■P AdS^ ^■■H 
a;" x^ x^ g 






/i, V, . . . 




indices 


g 






12 3 4 



The background geometry described so far is dual to ther- 
mal Af = 4 supersymmetric SU{Nc) Yang-Mills theory with 
Nf additional Af ~ 2 hypermultiplets. These hypermulti- 
plets arise from the lowest excitations of the strings stretching 
between the D7-branes and the background-generating D3- 
branes. The particles represented by the fundamental fields of 
the JV = 2 hypermultiplets model the quarks in our system. 
Their mass iriq is given by the asymptotic value of the sepa- 
ration of the D3- and D7-branes. In the coordinates used here 
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we write 



in 



2mq 

Vat 



M 

Y 



lim px{p) ~ 'Ti, 



(2.6) 



where we introduced the dimensionless scaled quark mass m. 

In addition to the parameters incorporated so far, we aim 
for a description of the system at finite chemical potential /i 
and baryon density ns- In field theory, a chemical poten- 
tial is given by a nondynamical time component of the gauge 
field. In the gravity dual, this is obtained by introducing a 
p-dependent gauge field component Aq{p) on the D7 brane 
probe. For now we consider a baryon chemical potential 
which is obtained from the U{1) subgroup of the flavor sym- 
metry group. The sum over flavors then yields a factor of Nf 
in front of the DBI action written down below. 

The value of the chemical potential p in the dual field the- 
ory is then given by 



I^L = lim Ao{p) 



Qh ~ 



2'Ka' 



(2.7) 



where we introduced the dimensionless quantity p for conve- 
nience. We apply the same normalization to the gauge field 
and distinguish the dimensionful quantity A from the dimen- 
sionless Aq = Aq {2TTa')/ Qh- 

The action for the probe branes' embedding function and 
gauge fields on the branes is 



5n 



(2.8) 



Here g is the induced metric ( 12.41 ) on the brane, F is the field 
strength tensor of the gauge fields on the brane and ^ are the 
branes' worldvolume coordinates. Tq? is the brane tension 
and the factor Nf arises from the trace over the generators of 
the symmetry group under consideration. For finite baryon 
density, this factor will be different from that at finite isospin 
density. 

In 181], the dynamics of this system of branes and gauge 
fields was analyzed in view of describing phase transitions at 
finite baryon density. Here we use these results as a starting 
point which gives the background configuration of the brane 
embedding and the gauge field values at finite baryon density. 
To examine vector meson spectra, we will then investigate the 
dynamics of fluctuations in this gauge field background. 

In the coordinates introduced above, the action S'dbi for the 
embedding x{p) ™^ the gauge fields' field strength F is ob- 
tained by inserting the induced metric and the field strength 
tensor into ( I2.8l l. As in iit], we get 



5dbi = -NfT^^Qlj /d«e ^//(l - x') 



X l/l-x'+p'x''- 2-^(1 -x')i^p'o , (2.9) 



where Fp^ = dpA^ is the field strength on the brane. 
depends solely on p. 



According to ||8|], the equations of motion for the back- 
ground fields are obtained after Legendre transforming the 
action (12.91 ). Varying this Legendre transformed action with 
respect to the field x gives the equation of motion for the em- 
beddings x(p)' 



- x")x^ L Sd^ 



2, /2 



3{l-x') + 2p'x 



2Ad 



2 1 - x' + p'x'' 



(2.10) 

The dimensionless quantity d is a constant of motion. It is 
related to the baryon number density jib by 



(2.11) 



Below, equation ( 12.10b will be solved numerically for differ- 
ent initial values xo d. The boundary conditions used are 



X(P = 1) = Xo, dpx{p) 



0. 



p=i 



(2.12) 



The quark mass m is determined by xo- It is zero for Xo = 
and tends to infinity for xo ^ 1. Figure [2] shows the de- 
pendence of the scaled quark mass m ~ 2mql\f\T on the 
starting value xo for different values of the baryon density 
parametrized by d cx 713. In general, a small (large) xo is 
equivalent to a small (large) quark mass. For xo < 0.5, xo 
can be viewed as being proportional to the large quark masses. 
At larger xo for vanishing d = 0, the quark mass reaches a fi- 
nite value. In contrast, at finite baryon density, if xo is close 
to 1, the mass rapidly increases when increasing xo further 
In embeddings with a phase transition, there exist more than 
one embedding for one specific mass value. In a small regime 
close to Xo = 1- there are more than one possible value of xo 
for a given m. So in this small region, xo is not proportional 

to TTlq. 

The equation of motion for the background gauge field A is 



dpA^ = 2d 



fVi~x' + p'x''' 



m'x')[p'p{i'x'r + 



(2.13) 



Integrating both sides of the equation of motion from pe to 
some p, and respecting the boundary condition Ao{p = 1) = 
lUt], we obtain the full background gauge field 

i ^f{i-x')[p'p{i-x^r+8d^] 

(2.14) 

Recall that the chemical potential of the field theory is given 
by limp^oo Aq{p) and thus can be obtained from the formula 
above. Examples for the functional behavior of Aq{p) are 
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FIG. 2: The dependence of the scaled quark mass m = 
on the horizon value xo = linip^i x of the embedding. 



shown in figure[3] Note that at a given baryon density 7^ 
there exists a minimal chemical potential which is reached in 
the limit of massles quarks. 

The asymptotic form of the fields x{p) Ao(/3) can be 
found from the equations of motion in the boundary limit p 

00, 



A, 







d Qh 



m c 
X = — ^ ^ 



(2.15) 
(2.16) 



Here p is the chemical potential, m is the dimensionless quark 
mass parameter given in ( 12.61 ), c is related to the quark conden- 
sate (but irrelevant in this work) and d is related to the baryon 
number density as stated in ( I2.11l i. See also figure[3]for this 
asymptotic behavior. The p-coordinate runs from the horizon 
value p = 1 to the boundary at p = 00. In most of this range, 
the gauge field is almost constant and reaches its asymptotic 
value, the chemical potential p, at p ^ 00. Only near the 
horizon the field drops rapidly to zero. For small xo ^ 0, 
the curves asymptote to the lowest (red) curve. So there is 
a minimal chemical potential for fixed baryon density in this 
setup. At small baryon density (d <C 0.00315) the embed- 
dings resemble the Minkowski and black hole embeddings 
known from the case without a chemical potential. Only a 
thin spike always reaches down to the horizon. 

In the setup described in this section we restrict ourselves 
to the regime of so called 'black hole embeddings' which are 
those embeddings ending on the horizon of the black hole, op- 
posed to 'Minkowski embeddings' , which would reach p = 
without touching the horizon. The black hole embeddings we 
use for this work (see figure O are not capable of describing 
matter in all possible phases. In fact we are able to cover the 
regime of fixed > and thus examine thermal systems in 
the canonical ensemble at finite baryon density. For a detailed 
discussion of this aspect see jgl fioll . 
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FIG. 3: The three figures of the left column show the embedding 
function x versus the radial coordinate p, the corresponding back- 
ground gauge fields Aq and the distance L — px between the D3 and 
the D7-branes at d = 10^^/4. L is plotted versus r, given by p^ = 
+ L^. In the right column, the same three quantities are depicted 
for d = 0.25. The five curves in each plot correspond to parametriza- 
tions of the quark mass to temperature ratio with xo ~ x(l) = 
0, 0.5, 0.9, 0.99 (all solid) and 0.99998 (dashed) from bottom up. 
These correspond to scaled quark masses m = 
0, 0.8089, 1.2886, 1.3030, 1.5943 in the left plot and to m = 
0, 0.8342, 1.8614, 4.5365, 36.4028 on the right. The curves on 
the left exhibit p « 10~*. Only the upper most curve on the 
left at Xo = 0.99998 develops a large chemical potential of p = 
0.107049. In the right column curves correspond to chemical po- 
tential values p = 0.1241, 0.1606, 0.5261, 2.2473, 25.3810 from 
bottom up. 



B. Holographic spectral functions 

Spectral functions contain information about the quasipar- 
ticle spectrum of a given theory. Recently, methods were de- 
veloped to compute spectral functions from the holographic 
duals of strongly coupled finite temperature gauge theories. 
In this work we extend these results to investigate the quasi- 
particle spectrum coiTesponding to vector mesons in the limit 
of vanishing spatial momentum. Therefore, we analyze the 
holographic dual to spectral functions for thermal J\f = 4 
supersymmetric SU{Nc) Yang-Mills theory with Nj funda- 
mental degrees of freedom (quarks) at finite baryon density 
and finite chemical potential. We compute the spectral densi- 
ties for the flavor current J, which is dual to the fluctuations 
A of the flavor gauge field on the supergravity side. 

Within field theory, the spectral function q) of some 
operator J{x) is defined via the imaginary part of the retarded 
Green function as follows 



(2.17) 
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where Energy w and spatial momentum q may be written in 
a four vector k = (w, q) and the Green function may be 
written as 

G'^{uj,q)^-i fd^xe'^^0{x°){[Jix),J{O)]) (2.18) 



One may find singularities of G^{uj, q) in the lower half of 
the complex cj-plane, including hydrodynamic poles of the 
retarded real-time Green function. Consider for example 



G' 



1 



These poles emerge as peaks in the spectral densities. 



5n = 



2r 



(2.19) 



(2.20) 



located at loq with a width given by F. These peaks are inter- 
preted as quasi-particles if their lifetime 1 /F is considerably 
long, i.e. if F ^ tJo- 

In this paper we use the gauge/gravity duality prescription 
of ifTill for calculating Green functions in Minkowski space- 
time. For further reference, we outline this prescription briefly 
in the subsequent. Starting out from a classical supergravity 
action 5ci for the gauge field A, according to |[l4ll we extract 
the function B{p) (containing metric factors and the metric 
determinant) in front of the kinetic term (9pA)^, 



(2.21) 



5ci= dpA^x B{p){dpAY 



Then we perform a Fourier transformation and solve the lin- 
earized equations of motion for the fields A in momentum 
space. The solutions in general are functions of all five co- 
ordinates in Anti-de Sitter space. Near the boundary we may 
separate the radial behavior from the boundary dynamics by 
writing 



A{p,k) = f{p,k)A'"'y{k). 



(2.22) 



where A'"'^(fc) is the value of the supergravity field at the 
boundary of AdS depending only on the four flat boundary 
coordinates. Thus by definition we have f{p, fc)|p_>tx) = 1- 
Then the retarded thermal Green function is given by 



G\uj,q) = 2B{p)J{p,~k)dpf{p,k) 



(2.23) 



The thermal correlators obtained in this way display hydro- 
dynamic properties, such as poles located at complex frequen- 
cies. They are used to compute the spectral densities ( 12.17b . 
We are going to compute the functions A{p, k) numerically in 
the limit of vanishing spatial momentum q ^ 0. The func- 
tions /(p, k) are then obtained by dividing out the boundary 
value A'"^^(fc) ~ limp^oo A{p, k). Numerically we obtain the 
boundary value by computing the solution at a fixed large p. 



III. SPECTRAL FUNCTIONS AT FINITE BARYON 
DENSITY 

A. Baryon diffusion 

In this section we calculate the baryon diffusion coefficient 
and its dependence on the baryon density. As discussed in 
ifioll . the baryon density affects the location and the pres- 
ence of the fundamental phase transition between two black 
hole embeddings observed in lUt]. This first order transition is 
present only very close to the separation line between the re- 
gions of zero and non-zero baryon density shown in figure [1] 

We show that this fundamental phase transition may also 
be seen in the diffusion coefficient for quark diffusion. In 
order to compute the diffusion using holography, we use the 
membrane paradigm approach developed in 12411 and extended 
in [17]. This method allows to compute various transport 
coefficients in Dp/Dq-brane setups from the metric coeffi- 
cients. The resulting formula for our background is the same 
as in d, 



D 



.911 V ^.900.944 



p=i 



dp 



'.900g44 



(3.1) 



The dependence of D on the baryon density and on the 
quark mass originates from the dependence of the embedding 
X on these variables. The results for D are shown in figure|4] 
The thick solid line shows the diffusion constant at vanishing 
baryon density found in in iflTll . which reaches D = at the 
fundamental phase transition. Increasing the baryon density, 
the diffusion coefficient curve is lifted up for small tempera- 
tures, still showing a phase transition up to the critical density 
d* = 0.00315. This is the same value as found in IS] in the 
context of the phase transition of the quark condensate. 

The diffusion coefficient never vanishes for finite density. 
Both in the Umit of T/M and T/M oo, D ■ T con- 
verges to l/27r for all densities, i.e. to the same value as for 
vanishing baryon density, as given for instance in ll24ll for R- 
charge diffusion. At the phase transition, the diffusion con- 
stant develops a nonzero minimum at finite baryon density. 
Furthermore, the location of the first order phase transition 
moves to lower values of T/M while we increase d towards 
its critical value. 

In order to give a physical explanation for this behavior, 
we focus on the case without baryon density first. We see 
that the diffusion coefficient vanishes at the temperature of the 
fundamental deconfinement transition. This is simply due to 
the fact that at and below this temperature, all charge carriers 
are bound into mesons not carrying any baryon number. 

For non-zero baryon density however, there is a fixed num- 
ber of charge carriers (free quarks) present at any finite tem- 
perature. This implies that the diffusion coefficient never van- 
ishes. Switching on a very small baryon density, even be- 
low the phase transition, where most of the quarks are bound 
into mesons, by definition there will still be a finite amount 
of free quarks. By increasing the baryon density, we increase 
the amount free quarks, which at some point outnumber the 
quarks bound in mesons. Therefore in the large density limit 
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0.775 



FIG. 4: The diffusion coefficient times temperature is plotted 
against the mass-scaled temperature for diverse baryon densities 
parametrized by d — 0.1 (uppermost line in upper plot, not visible in 
lower plot), 0.004, (long-dashed), 0.00315 (thin solid), 0.002 (long- 
short-dashed), 0.000025 (short-dashed) and (thick solid). The fi- 
nite baryon density lifts the curves at small temperatures. Therefore 
the diffusion constant never vanishes but is only minimized near the 
phase transition. The lower plot zooms into the region of the tran- 
sition. The phase transition vanishes above a critical value d* = 
0.00315. The position of the transition shifts to smaller T/M, as d 
is increased towards its critical value. 



the diffusion coefficient approaches = 1/ {2ttT) for all 
values of T/M, because only a negligible fraction of the 
quarks is still bound in this limit. 

Note that as discussed in IH |^ [l^ there exists a region 
in the {nB,T) phase diagram at small ub and T where the 
embeddings are unstable. In figure]?] this corresponds to the 
region just below the phase transition at small baryon density. 
This instability disappears for large ns- 



B. Vector mesons in the black hole phase 

1. Application of calculation method 

We now compute the spectral functions of flavor currents 
at finite baryon density ns, chemical potential /i and tem- 
perature in the 'black hole phase'. As black hole phase the 
authors of denote the phase of matter which has nonzero 
baryon density. Cornpared to the limit of vanishing chemical 
potential treated in iflTll . we discover a qualitatively different 
behavior of the finite temperature oscillations corresponding 
to vector meson resonances. 



To obtain the spectral functions, we compute the correla- 
tions of flavor gauge field fluctuations about the back- 
ground given by i2.9i . denoting the full gauge field by 

= S^Mp) + A^{x,p) . (3.2) 



According to section III Al the background field has a non- 
vanishing time component, which depends solely on p. The 
fluctuations in turn are gauged to have non-vanishing compo- 
nents along the Minkowski coordinates x only and only de- 
pend on these coordinates and on p. Additionally they are 
assumed to be small, so that it suffices to consider their lin- 
earized equations of motion. At this point we neglect the fluc- 
tuation of the scalar and pseudoscalar modes and their cou- 
pling to the vector fluctuations. In fact there is no such cou- 
pling in the limit of vanishing spatial momentum, which we 
restrict to below. 

The resulting equations of motion are obtained from the ac- 
tion (I2.8I 1, where we introduce small fluctuations A by setting 
^ = 2 A^] with A ^ A + A. The background 
gauge field A is given by ( 12.131 ). The fluctuations now propa- 
gate on a background G given by 



G = g + F, 

and their dynamics is determined by the Lagrangian 



£= v/|det(G + F)|, 



(3.3) 



(3.4) 



with the fluctuation field strength = 2d^^A^y Since the 
fluctuations and their derivatives are chosen to be small, we 
consider their equations of motion only up to linear order, as 
derived from the part of the Lagrangian C which is quadratic 
in the fields and their derivatives. Denoting this part by C2, 
we get 



1 



£2 =--V|detG| 



(3.5) 



Here and below we use upper indices on G to denote elements 
of G^^. The equations of motion for the components of A are 



0=9, 



VidetGI 



X ( G^"'G''^ - G^"'G''''-G^'"'^G^^'] d^^A^,^ 



(3.6) 



The terms of the corresponding on-shell action at the p- 
boundaries are (with p as an index for the coordinate p, not 
summed) 



00 

On 



( {G°^fAodpAo - G^^G'^AdpAk 



AoG^"tr(G-^F) 



(3.7) 
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Note that on the boundary ps p ^ oo, the background ma- 
trix G reduces to the induced D7-brane metric g. Therefore, 
the analytic expression for boundary contributions to the on- 
shell action is identical to the one found in iITtIi . There, the co- 
ordinates in Minkowski directions were chosen such that the 
fluctuation four vector k exhibits only one non vanishing spa- 
tial component, e. g. in x-direction as A: = (oj, 0, 0). Then 
the action was expressed in terms of the gauge invariant field 
component combinations 



Err = UjAr. 



E,, 



(3.8) 



In the case of vanishing spatial momentum g ^ 0, the Green 
functions for the different components coincide and were 
computed as iIitIi 



Gi 



Gi 



NfNcT^ 
8 



lim 

p — >oo 



E{p) J ' 
(3.9) 

where the E{p) in the denominator divides out the boundary 
value of the field in the limit of large p, as discussed after 
( I2.23l l. The indices on the Green function denote the compo- 
nents of the operators in the correlation function, in our case 
all off-diagonal correlations (as Gyz, for example) vanish. 

In our case of finite baryon density, new features arise 
through the modified embedding and gauge field background, 
which enter the equations of motion ( I3.6l l for the field fluctu- 
ations. To apply the prescription to calculate the Green func- 
tion, we Fourier transform the fields as 



e^''^A^{p,k). 



(3.10) 



As above, we are free to choose our coordinate system to 
give us a momentum vector of the fluctuation with nonvanish- 
ing spatial momentum only in .x-direction, k = {oj, q, 0, 0). 

For simplicity we restrict ourselves to vanishing spatial mo- 
mentum q = 0. In this case the equations of motion for 
transversal fluctuations Ey^z match those for longitudinal fluc- 
tuations Ex- For a more detailed discussion see |fl7||. As an 
example consider the equation of motion obtained from (13.6b 
with CT = 2, determining Ey = L0A2, 



0=E 



„ , d,[^\d^\G'^G^^] G-oo 2 2p 



v/|dctG|G22G44 



-E" + 



dp In 



■E 




The symbol tr denotes the dimensionless frequency tr ~ 
ijj/{2ttT), and we made use of the dimensionless radial co- 
ordinate p. 

In order to numerically integrate this equation, we deter- 
mine local solutions of that equation near the horizon p = 1 . 



These can be used to compute initial values in order to inte- 
grate ( 13.111 ) forward towards the boundary. The equation of 
motion ( 13.1 II ) has coefficients which are singular at the hori- 
zon. According to standard methods jzsll . the local solution 
of this equation behaves as [p — pu)^, where /3 is a so-called 
'index' of the differential equation. We compute the possible 
indices to be 



l3 = ±ito. 



(3.12) 



Only the negative one will be retained in the following, since 
it casts the solutions into the physically relevant incoming 
waves at the horizon and therefore satisfies the incoming wave 
boundary condition. The solution E can be split into two fac- 
tors, which are {p — l)"'"" and some function F{p), which is 
regular at the horizon. The first coefficients of a series expan- 
sion of can be found recursively as described in lfl5HT6ll . 
At the horizon the local solution then reads 



Eip)^ip-1)- 



F{P) 



(3.13) 



So, F{p) asymptotically assumes values 



F{p=l) = l, 



dpF{p) 



"2 



= -r- (3.14) 



For the calculation of numbers, we have to specify the 
baryon density d and the mass parameter xo ^ 'mq/T to ob- 
tain the embeddings x used in ( 13.1 lb . Then we obtain a so- 
lution for a given frequency ro using initial values ( 13. 13b and 
(13.14b in the equation of motion ( 13.111 ). This eventually gives 
us the numerical solutions for E{p). 

Spectral functions are then obtained by combining ( 13.91 ) and 
dTlTI i, 



9^(w,0) 



NfNcT^ 



Im lim 

p — ^00 



. dpEjp) 
E{p) 



(3.15) 



2. Results for spectral functions 

We now discuss the resulting spectral functions at finite 
baryon density, and observe crucial qualitative differences 
compared to the case of vanishing baryon density. In fig- 
ures |5] to [8] some examples for the spectral function at fixed 
baryon density ub oc d are shown. To emphasize the reso- 
nance peaks, in some plots we subtract the quantity 



(3.16) 



around which the spectral functions oscillate, cf. figure |9l 

The graphs are obtained for a value of d above d* (given by 
dl.lb ). where the fundamental phase transition does not occur. 
The different curves in these plots show the spectral functions 
for different quark masses, corresponding to different posi- 
tions on the solid blue line in the phase diagram shown in fig- 
ure[T] Regardless whether we chose d to be below or above the 
critical value of d, we observe the following behavior of the 
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FIG. 5: The finite temperature part of the spectral function $H— 5Ho (in 
units of NfNcT'^/4) at finite baryon density d. The maximum grows 
and shifts to smaller frequencies as xo is increased towards xo = 
0.7, but then turns around to approach larger frequency values. 
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FIG. 6: The finite temperature part of the spectral function — ?Ho 
(in units of Nf NcT^ /^) at finite baryon density d. In the regime 
of Xo shown here, the peak shifts to larger frequency values with 
increasing xo- 



spectral functions with respect to changes in the quark mass 
to temperature ratio. 

Increasing the quark mass from zero to small finite values 
results in more and more pronounced peaks of the spectral 
functions. This eventually leads to the formation of resonance 
peaks in the spectrum. At small masses, though, there are no 
narrow peaks. Only some maxima in the spectral functions 
are visible. At the same time as these maxima evolve into 
resonances with increasing quark mass, their position changes 
and moves to lower freqencies in, see figure|5] This behavior 
was also observed for the case of vanishing baryon density in 
d. 

However, further increasing the quark mass leads to a cru- 
cial difference to the case of vanishing baryon density. Above 
a value m'"'" of the quark mass, parametrized by Xg™, the 
peaks change their direction of motion and move to larger val- 
ues of to, see figure|6l Still the maxima evolve into more and 
more distinct peaks. 
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FIG. 7: The finite temperature part ?H — SHo of the spectral func- 
tion (in units of NfNcT^/^) at finite baryon density d. The oscilla- 
tion peaks narrow and get more pronounced compared to smaller xo. 
Dashed vertical lines show the meson mass spectrum given by equa- 
tion frn\ . 



Eventually at very large quark masses, given by x closer 
and closer to 1, the positions of the peaks asymptotically reach 
exactly those frequencies which correspond to the masses of 
the vector mesons at zero temperature lfl9l] . In our coordi- 
nates, these masses are given by 



M 



-^V2(« + l)(n + 2). 



(3.17) 



where labels the Kaluza-Klein modes arising from the D7- 
brane wrapping S^, and ioo is the radial distance in (8,9)- 
direction between the stack of D3-branes and the D7, evalu- 
ated at the AdS'-boundary, 



lim gx{g)- 



(3.18) 



The formation of a line-like spectrum can be interpreted as 
the evolution of highly unstable quasi-particle excitations in 
the plasma into quark bound states, finally turning into nearly 
stable vector mesons, cf. figures|7]and[8] 

We now consider the turning behavior of the resonance 
peaks shown in figures |5] and |6] There are two different sce- 
narios, depending on whether the quark mass is small or large. 

First, when the quark mass is very small niq <C T, we are 
in the regime of the phase diagram corresponding to the right 
half of figure[T] In this regime the influence of the Minkowski 
phase is negligible, as we are deeply inside the black hole 
phase. We therefore observe only broad structures in the spec- 
tral functions, instead of peaks. 

Second, when the quark mass is very large, rUg ^ T, or 
equivalently the temperature is very small, the quarks behave 
just as they would at zero temperature, forming a line-like 
spectrum. This regime corresponds to the left side of the phase 
diagram in figure [T] where all curves of constant d asymptote 
to the Minkowski phase. 

The turning of the resonance peaks is associated to being in 
the first or in the second regime. At Xq™ the two regimes are 
connected and none of them is dominant. 
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FIG. 8: The spectral function $R (in units of MfNcT'^/A) at finite 
baryon density d. At large xo, as here, the peaks approach the dashed 
drawn line spectrum given by ( 13.17b . 



60000 



d = 0.25 




50000- Xo= 0.997 

: fR(lT),0) 

40000 - 

IHo 

30000 : 
20000 
10000 




FIG. 9: The thermal spectral function (in units of NfNcT'^/4) 
compared to the zero temperature result IHo. 



The turning behavior is best understood by following a line 
of constant density d in the phase diagram of figure [1] Con- 
sider for instance the solid blue line in figure [T] starting at 
large temperatures/small masses on the right of the plot. First, 
we are deep in the unshaded region {ns ^ 0), far inside the 
black hole phase. Moving along to lower T/ AI, the solid blue 
line in figure[T]rapidly bends upwards, and asymptotes to both 
the line corresponding to the onset of the fundamental phase 
transition, as well as to the separation line between black hole 
and Minkowski phase (gray region). 

This may be interpreted as the quarks joining in bound 
states. Increasing the mass further, quarks form almost sta- 
ble mesons, which give rise to resonance peaks at larger fre- 
quency if the quark mass is increased. The confined and de- 
confined phase are coexistent asymptotically for T/ M 0. 

We also observe a dependence of Xo"" on the baryon den- 
sity. As the baryon density is increased from zero, the value 
of Xo™ decreases. 

Figures [8] and |9] show that higher n excitations from the 
Kaluza-Klein tower are less stable. While the first resonance 
peaks in this plot are very narrow, the following peaks show a 
broadening with decreasing amplitude. 

This broadening of the resonances is due to the behaviour 
of the quasinormal modes of the fluctuations, which corre- 
spond to the poles of the correlators in the complex u; plane, 
as described in the example ( |2.19l l and sketched in figure [TO] 
The location of the resonance peaks on the real frequency axis 
corresponds to the real part of the quasinormal modes. It is a 
known fact that the the quasinormal modes develop a larger 
real and imaginary part at higher n. So the sharp resonances 
at low ro, which correspond to quasiparticles of long lifetime, 
originate from poles whith small imaginary part. For higher 
excitations in n at larger ro, the resonances broaden and get 
damped due to larger imaginary parts of the corresponding 
quasi normal modes. 

For increasing mass we described above that the peaks of 
the spectral functions first move to smaller frequencies un- 
til they reach the turning point m'"'". Further increasing the 
mass leads to the peaks moving to larger frequencies, asymp- 




FIG. 10: Qualitative relation between the location of the poles in 
the complex frequency plane and the shape of the spectral functions 
on the real uj axis. The function plotted here is an example for the 
imaginary part of a correlator. Its value on the real u axis represents 
the spectral function. The poles in the right plot are closer to the real 
axis and therefore there is more structure in the spectral function. 



totically approaching the line spectrum. This behavior can be 
translated into a movement of the quasinormal modes in the 
complex plane. It would be interesting to compare our results 
to a direct calculation of the quasinormal modes of vector fluc- 
tuations in analogy to i^. 

In [6] the quasinormal modes are considered for scalar fluc- 
tuations exclusively, at vanishing baryon density. The au- 
thors observe that starting from the massless case, the real 
part of the quasinormal frequencies increases with the quark 
mass first, and then turns around to decrease. This behavior 
agrees with the peak movement for scalar spectral functions 
observed in IitI figure 9] (above the fundamental phase tran- 
sition, Xo < 0.94) where the scalar meson resonances move 
to higher frequency first, turn around and move to smaller fre- 
quency increasing the mass further. These results do not con- 
tradict the present work since we consider vector modes ex- 
clusively. The vector meson spectra considered in UtIi at van- 
ishing baryon density only show peaks moving to smaller fre- 
quency as the quark mass is increased. Note that the authors 
there continue to consider black hole embeddings below the 
fundamental phase transition which are only metastable, the 
Minkowski embeddings being thermodynamic ally favored. 
At small baryon density and small quark mass our spectra 
are virtually coincident with those of ITtIi . In our case, at 
finite baryon density, black hole embeddings are favored for 
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all values of the mass over temperature ratio. At small val- 
ues of T/M in the phase diagram of figure [1] we are very 
close to the Minkowski regime, temperature effects are small, 
and the meson mass is proportional to the quark mass as in 
the supersymmetric case. Therefore, the peaks in the spectral 
function move to the right (higher frequencies) as function of 
increasing quark mass. 

The turning point in the location of the peaks is a con- 
sequence of the transition between two regimes, i.e. the 
temperature-dominated one also observed in IitIi . and the 
potential-dominated one which asymptotes to the supersym- 
metric spectrum. 

We expect the physical interpretation of the left-moving of 
the peaks in the temperature-dominated regime to be related 
to the strong dissipative effects present in this case. This is 
consistent with the large baryon diffusion coefficient present 
in this regime as discussed in section 1111 Al and shown in fig- 
ure|4] A detailed understanding of the physical picture in this 
regime requires a quantitative study of the quasipaticle behav- 
ior which we leave to future work. 

In our approach it is straightforward to investigate the T 
limit since black hole embeddings are thermodynamically 
favored even near T = at finite baryon density. We expect 
that a right-moving of the peaks consistent with the SUSY 
spectrum should also be observable for Minkowski embed- 
dings at vanishing baryon density. However this has not been 
investigated for vector modes neither in ^ nor in fl7^. 



IV. SPECTRAL FUNCTIONS AT FINITE ISOSPIN 
DENSITY 

A. Radially varying SU (2)-background gauge field 



field strength, such that we get two copies of the baryonic 
background Fpo on the diagonal of the flavor matrix. 



dpAo 






'dpAo 



(4.4) 



The action for the isospin background differs from the ac- 
tion i2.9i for the baryonic background only by a group the- 
oretical factor: The factor Tr = 1/2 (compare ( 14. It ) replaces 
the baryonic factor Nf in equation ( 12. 8t , which arises by sum- 
mation over the U{1) representations. We can thus use the 
embeddings x{p) ™d background field solutions Ao{p) of 
the baryonic case of fsfl, listed here in section UTAI As be- 
fore, we collect the induced metric g and the background field 
strength F in the background tensor G ^ g + F. 

We apply the background field method in analogy to the 
baryonic case examined in section |III] As before, we obtain 
the quadratic action by expanding the determinant and square 
root in fluctuations A^. The term linear in fluctuations again 
vanishes by the equation of motion for our background field. 
This leaves the quadratic action 



oo 

= QH{27T^R')TD7Tr Jdpd^X ^/\d^\ 



1 

4 



QH 



(27ra 



Qh 



(4.5) 



In order to examine the case Nf = 2 in the strongly cou- 
pled plasma, we extend our previous analysis of vector meson 
spectral functions to a chemical potential with SU (2)-flavor 
(isospin) structure. Starting from the general action 



T.T^T /d«e Vldet(ff + i^)|, (4.1) 



we now consider field strength tensors 



Qh j^abc 2b AC \ 

2W^ ^^^^ ) ' 



(4.2) 



with the Pauli matrices cr" and A given by equation (13. 2t . The 
factor £i|^/(27ra') is due to the introduction of dimensionless 
fields as described below (12.7b . In order to obtain a finite 
isospin-charge density nj and its conjugate chemical poten- 
tial pi, we introduce an SU (2)-background gauge field A ll22ll 



Ay 





-1 



(4.3) 



This specific choice of the 3-direction in flavor space as well 
as spacetime dependence simplifies the isospin background 



Note that besides the familiar Maxwell term, two other terms 
appear, which are due to the non-Abelian structure. One of 
the new terms depends linearly, the other quadratically on the 
background gauge field A and both contribute nontrivially to 
the dynamics. The equation of motion for gauge field fluctua- 
tions on the D7-brane is 



Vld^^ {a'^G"" - G'^'G'"') f;^^ 



(4.6) 



V|detG||^igr''3 {G''°G'"' - G'""G°'^) 



27ra 



with the modified field strength linear in fluctuations = 
2a[^^^] + r''Al{5o,At + <5o.4)^.|,/(2W). 

Integration by parts of (14.5b and application of ( 14.6b yields 
the on-shell action 



on-shell 
iso 



QHTrTDTTT^R^ Jd*x VldetGI 



Qi'iQu't^ _ G-'^'G^^n Al,F'^^ 



Pb 

PH 



(4.7) 



The three flavor field equations of motion (flavor index a 
1,2,3) for fluctuations in transversal Lorentz-directions a 
2, 3 can again be written in terms of the combination i?^ 
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qAq + ujA'^. At vanishing spatial momentum 5 = we get 

(4.8) 







V|detG|G44G22 ^ 



V27ra' 



27ra 



= 



2// , dpi^\dctG\G^^G^' 
v/|detG|G44G22 



(4.9) 



G°" 

G44 



■ ^0 



V27ra' 



2zeHC^G00 el, .3 



G44 27ra' 



= 



9p(^d^G^4G'- 
\/|detG|G44G22 



-E^ — 



G44 



(4.10) 



Note that we use the dimensionless background gauge 
field Al = Al{2TTa')/gH and qh = ttTR^. Despite the 
presence of the new non-Abelian terms, at vanishing spatial 
momentum the equations of motion for longitudinal fluctua- 
tions are the same as the transversal equations ( 14.8b . ( 14. 9t and 
(|4TT0] |. such that E ^ Et = El. 

Note at this point that there are two essential differences 
which distinguish this setup from the approach with a constant 
potential A^ at vanishing mass followed in ll22ll . First, the in- 
verse metric coefficients g*^^ contain the embedding function 
x{p) computed with varying background gauge field. Second, 
the background gauge field Ag giving rise to the chemical po- 
tential now depends on p. 

Two of the ordinary second order differential equa- 
tions ( |4.8t . ( 14.91 ). ( 14.101 ) are coupled through their flavor sti'uc- 
ture. Decoupling can be achieved by transformation to the 
flavor combinations |22l 



X = E^ + iE^, Y ^ E^ -iE'^ . 
The equations of motion for these fields are given by 



^|detG|G44G22 



(4.11) 



(4.12) 



Q_y. , ap(^d^G^^G^^) ^^, 
v/|detG|G44G22 



(4.13) 



( IT) 

\R4 G44 ^"^ 



xnyY. 



= E- 



,/|detG|G44G22 



i?4 G44 



(4.14) 



with dimensionless m = Al/{2TrT) and ro = w/(27rT). Pro- 
ceeding as described in section|lIll we determine the local so- 
lution of (I4.I2I 1. (14.131 1 and (14.14b at the horizon. The indices 
turn out to be 



(27rT) 



(4.15) 



Since Aq{p = 1) = in the setup considered here, we are left 
with the same index as in (13.12b for the baryon case. There- 
fore, here the chemical potential does not influence the sin- 
gular behavior of the fluctuations at the horizon. The local 
solution coincides to linear order with the baryonic solution 
given in ( 13.131 ). 

Application of the recipe described in section IllBI yields 
the spectral functions of flavor current correlators shown in 
figures [TT] and [12] Note that after transforming to flavor com- 
binations X and Y, given in ( 14.1 lb . the diagonal elements 
of the propagation submatrix in flavor-transverse X, Y direc- 
tions vanish, Gxx = Gyy — 0, while the off-diagonal el- 
ements give non- vanishing contributions. However, the com- 
ponent E^, longitudinal in flavor space, is not influenced by 
the isospin chemical potential, such that G^js is nonzero, 
while other combinations with E'^ vanish (see 1I22I1 for de- 
tails). 

Introducing the chemical potential as described above for 
a zero-temperature AdS^ x background, we obtain the 
gauge field correlators in analogy to ll26ll . The resulting spec- 
tral function for the field theory at zero temperature but finite 
chemical potential and density £Ho,iso is given by 



NrT^Tr 



47r(tr ± moo)^ 



(4.16) 



with the dimensionless chemical potential rrioo = 
limp^oo Ay{2TTT) = /i/(27rT). Note that ( l4l6l l is in- 
dependent of the temperature. This part is always subtracted 
when we consider spectral functions at finite temperature, in 
order to determine the effect of finite temperature separately, 
as we did in the baryonic case. 



B. Results at finite isospin density 

In figure [TT] we compare typical spectral functions found 
for the isospin case (solid lines) with that found in the bary- 
onic case (dashed line). While the qualitative behavior of the 
isospin spectral functions agrees with the one of the bary- 
onic spectral functions, there nevertheless is a quantitative dif- 
ference for the components X, Y, which are transversal to 
the background in flavor space. We find that the propagator 
for flavor combinations Gyx exhibits a spectral function for 
which the zeroes as well as the peaks are shifted to higher fre- 
quencies, compared to the Abelian case curve. For the spec- 
tral function computed from Gxy, the opposite is true. Its 
zeroes and peaks appear at lower frequencies. As seen from 
figure [12] also the quasiparticle resonances of these two dif- 
ferent flavor correlations show distinct behavior The quasi- 
particle resonance peak in the spectral function D\yx appears 
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FIG. 1 1 ; The finite temperature part of spectral functions fHiso — 
SHo,iso (in units of NcT^Tr/4:) of currents dual to fields X, Y are 
shown versus to. The dashed line shows the baryonic chemical po- 
tential case, the solid curves show the spectral functions in presence 
of an isospin chemical potential. Plots are generated for xo ~ 0.5 
and d = 0.25. The combinations XY and YX split in opposite 
directions from the baryonic spectral function. 
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FIG. 12: A comparison between the finite temperature part of the 
spectral functions ^ixY and SHi x (solid lines) in the two flavor di- 
rections X and Y transversal to the chemical potential is shown in 
units of NcT'^Tr/4: for large quark mass to temperature ratio xo = 
0.99 and d = 0.25. The spectral function iH^s^s along the a — 3- 
flavor direction is shown as a dashed line. We observe a splitting of 
the line expected at the lowest meson mass at tt) = 4.5360 (n = 0). 
The resonance is shifted to lower frequencies for SHxy and to higher 
ones for diyx, while it remains in place for JH^js^a- The second 
meson resonance peak (n — 1) shows a similar behavior So the 
different flavor combinations propagate differently and have distinct 
quasiparticle resonances. 



at higher frequencies than expected from the vector meson 
mass formula il.2\ (shown as dashed grey vertical lines in fig- 
ure [T2]i. The other flavor-transversal spectral function SKxy 
displays a resonance at lower frequency than observed in the 
baryonic curve. The spectral function for the third flavor di- 
rection y^E's behaves as E in the baryonic case. 

This may be viewed as a splitting of the resonance peak into 
three distinct peaks with equal amplitudes. This is due to the 
fact that we explicitly break the symmetry in flavor space by 
our choice of the background field A^. Decreasing the chem- 
ical potential reduces the distance of the two outer resonance 
peaks from the one in the middle and therefore the splitting is 



reduced. 

The described behavior resembles the mass splitting of 
mesons in presence of a isospin chemical potential expected to 
occur in QCD |!27','28'l. A linear dependence of the separation 
of the peaks on the chemical potential is expected. Our ob- 
servations confirm this behavior. Since our vector mesons are 
isospin triplets and we break the isospin symmetry explicitly, 
we see that in this respect our model is in qualitative agree- 
ment with effective QCD models. Note also the complemen- 
tary discussion of this point in ll29ll . 

To conclude this section, we comment on the relation of 
the present results to those of our previous paper i22ll where 
we considered a constant non-Abelian gauge field background 
for zero quark mass. From equation ( |4.15t . the difference be- 
tween a constant non-vanishing background gauge field and 
the varying one becomes clear In II22II the field is chosen 
to be constant in p and terms quadratic in the background 
gauge field Aq ^ 1 are neglected. This implies that the 
square {vo^mY in ( 14.12b and ( 14.13b is replaced by Xo'^^2K)m, 

such that we obtain the indices B = ±txiA/lT ."^w^. ' in- 



Stead of ( 14.15b . If we additionally assume to <C Aq, then the 1 
under the square root can be neglected. In this case the spec- 
tral function develops a non-analytic structure coming from 
the y/uj factor in the index. 

However in the case considered here, the background gauge 
field is a non-constant function of p which vanishes at the hori- 
zon. Therefore the indices have the usual form /3 = ±icj from 
(14. 15b . and there is no non-analytic behavior of the spectral 
functions, at least none originating from the indices. 

It will also be interesting to consider isospin diffusion in 
the setup of the present paper However, in order to see non- 
Abelian effects in the diffusion coefficient, we need to give 
the background gauge field a more general direction in fla- 
vor space or a dependence on further space-time coordinates 
besides p. In that case, we will have a non-Abelian term 
in the background field strength F^u = d^A'^ — d^A'^ + 
jabcj^b Ac^^^2 1 (^2Tia') in Contrast to dpAf^ considered here. 



V. CONCLUSION 

Two distinct setups were examined here at non-zero charge 
density in the black hole phase. First, switching on a baryon 
chemical potential at non-zero baryon density, we find that 
nearly stable vector mesons exist close to the transition line 
to the Minkowski phase. Far from this line, at small quark 
masses, we essentially recover the picture given in the case of 
vanishing chemical potential iITtIi . Increasing the quark mass 
beyond a distinct value, the plasma changes its behavior in or- 
der to asymptotically behave as it would at zero temperature. 
In the spectral functions we computed, this zero-temperature- 
like behavior is found in form of line-like resonances, ex- 
actly reproducing the zero-temperature supersymmetric vec- 
tor meson mass spectrum. A turning point at to = m'"™, 
where to = M /T, is observed: Below to'"™, the resonance 
peaks move to lower frequencies as function of rising quark 
mass. This behavior of the system resembles the behav- 
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ior known of that system without a chemical potential iHtIi . 
Above to'"™, the resonance peaks move to higher frequencies 
as function of the quark mass. This is the zero-temperature- 
like regime. Moreover, an examination of the diffusion coeffi- 
cient reveals that the phase transition separating two different 
black hole phases iQ] is shifted towards smaller temperature 
as the baryon density is increased. 

Second, we switched on a nonzero isospin density, and 
equivalently an isospin chemical potential arises. The spec- 
tral functions in this case show a qualitatively similar be- 
havior as those for baryonic potential. However, we ad- 
ditionally observe a splitting of the single resonance peak 
at vanishing isospin potential into three distinct resonances. 
This suggests that by explicitly breaking the flavor symme- 
try by a chemical potential, the isospin triplet states, vector 
mesons in our case, show a mass splitting similar to that ob- 
served for QCD 1270 . It is an interesting task to explore the 
features of this isospin theory in greater detail in order to 
compare with available lattice data and effective QCD mod- 
els H HI [H [M [M HI IM H?! HI . In most of these 
approaches, baryon and isospin chemical potential are con- 
sidered at the same time, which suggests another promising 
extension of this work. Moreover, in the context of gravity du- 
als, it will be interesting to compare our results for the isospin 
chemical potential to the recent work lEgll . 

Alternatively, instead of giving the gauge field time compo- 
nent a non-vanishing vev, one may also switch on B-field com- 
ponents and connect the framework developed in lf39[ l40l l4lll 
with the calculation of spectral functions for the dual gauge 
theory. 
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APPENDIX A: NOTATION 

The five-dimensional AdS Schwarzschild black hole space 
in which we work is endowed with a metric of signature 
( — ,+ ,+,+,+), as given explicitly in ( 12.11 ). We make use 
of the Einstein notation to indicate sums over Lorentz indices, 
and additionally simply sum over non-Lorentz indices, such 
as gauge group indices, whenever they occur twice in a term. 

To distinguish between vectors in different dimensions of 
the AdS space, we use bold symbols like q for vectors in the 
the three spatial dimensions which do not live along the radial 
AdS coordinate. Four-vectors which do not have components 
along the radial AdS coordinates are denoted by symbols with 
an arrow on top, as q. 

The Green functions G = {J I) considered give correla- 
tions between currents J and /. These currents couple to 
fields A and B respectively. In our notation we use symbols 
such as Gj^ajib to denote correlators of currents coupling to 

fields Al and Af, with flavor indices a, b and Lorentz indices 
k,l = 0, 1, 2, 3. If no other indices are of relevance for the 
discussion we restrict ourselves to Lorentz indices. For the 
gauge field combinations X and Y given in ( I4.11l i. we obtain 
Green functions Gxy or Gyx denoting correlators of the cor- 
responding currents. 
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